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ABSTRACT 

The advection of a passive scalar by incompressible turbulence is considered using recur- 
sive renormalization group procedures in the differential subgrid shell thickness limit. It is 
shown explicitly that the higher order nonlinearities induced by the recursive renormaliza- 
tion group procedure preserve Galilean invariance. Differential equations, valid for the entire 
resolvable wavenumber k range, are determined for the eddy viscosity and eddy diffusivity 
coefficients and it is shown that higher order nonlinearities do not contribute as k — ► 0, but 
have an essential role as k — > k c , the cutoff wavenumber separating the resolvable scales from 
the subgrid scales. The recursive renormalization transport coefficients and the associated 
eddy Prandtl number are in good agreement with the ^-dependent transport coefficients 
derived from closure theories and experiments. 
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1 Introduction 


The turbulent transport of a passive scalar, while serving sound pedagogical purposes, is 
also of interest in the spreading of temperature, humidity and pollution in the atmosphere 
as well as in other problems (Csanady 1 ). Here w r e shall apply recursive renormalization group 
(RNG) procedures to the subgrid modeling of a passive scalar field T(k,£) being advected 
by a turbulent Navier-Stokes velocity field u(k,£). Subgrid modeling is necessary for the 
high-Reynolds number turbulent flows of interest because of the limitations of current and 
foreseeable supercomputers 2 . Another advantage of considering the problem of passive scalar 
transport is that the spectral transport coefficients (eddy diffusivity and eddy viscosity) 
determined from our RNG theory can be compared to those arising from closure-based 
theories 3 ’ 4 . It should be noted that the transport coefficients in these closure theories are 
determined over the whole resolvable scales. 

Recently, two distinct approaches of RNG to fluid turbulence have arisen: one based on 
the work of Forster et. al. 5 , called e-RNG, and the other based on Rose 6 , called recursive- 
RNG. Some aspects of these two approaches have been discussed'. In particular, we point 
out here that in c-RNG, a small parameter t is introduced through the forcing correlation 
function. Yakhot & Orszag 8 have extrapolated e <C 1 to c — » 4 in order to reproduce 
the Kolmogorov energy spectrum. Furthermore, it is also necessary to take the distant 
interaction limit 9 , k — > 0. Thus, it is difficult to compare the transport coefficients generated 
by Kraichnan 3 and Chollet 4 , with that determined from e-RNG. 

In this paper, we continue our application of recursive RNG 10-11 to turbulence. The 
basic differences between the two RNG procedures are that in recursive RNG: 

(i) The e-expansion is not applied. 

(ii) The turbulent transport coefficients are determined for the whole resolvable wavenum- 
ber scales, 

(iii) Higher order nonlinearities are generated in the renormalized momentum equation. 

In Sec. 2 we derive the renormalized evolution equations for the passive scalar T(k, f) 
and the fluid velocity u(k, t ) as well as the recursion relations from which the eddy diffusivity 
and eddy viscosity can be determined. Because of the presence of higher order nonlinearities 
in the renormalized equations, it is not apparent that Galilean invariance is still preserved. 
These questions are addressed in Sec. 3 where we prove that the RNG-evolution equations are 
Galilean invariant: a property deemed necessary in any subgrid model 12 . It has been found 
very difficult to find fixed points for the RNG-difference recursion relations if the subgrid 
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shell thickness is chosen too small 1011 . If recursion RNG procedures are to be employed 
successfully in more complicated flow problems, then it is necessary that these difference 
recursion relations be simplified. In Sec. 4 we proceed to the differential limit of these 
recursion relations, paying careful attention to the k -> 0 limit. We show that the higher 
order RNG-induced nonlinearities do not contribute to the k — > 0 limit of the RNG recursion 
relations, but play a significant role for k — > k c , where k c is the wavenumber that separates the 
resolvable scale from the subgrid scale. In the Appendix, we contrast our ordinary differential 
equations for the RNG eddy viscosity and diffusivity with that generated by Yakhot-Orszag 8 
by their c-RNG, a theory that is valid only in the distant interaction limit k ► 0. In contrast, 
in the recursive RNG approach, the differential equations for the eddy transport coefficients 
are valid over the whole resolvable scales and not just at k ► 0, and no c expansion is needed. 
The turbulent transport coefficients for the second moments (i.e., for the time evolution of 
U a p(k,t) =< u a (k, t)up(—k,t) > and the scalar variance 0(k,Z) =< T(k, t)T(-k, t) > ) 
are determined in Sec. 5. In particular, the importance of the RNG-induced higher order 
nonlinearities is very apparent. The spectral eddy viscosity, diffusivity and Prandtl number 
are derived in Sec. 6, while in Sec. 7 we present our conclusions. 

2 Renormalized momentum equation for velocity and 
passive scalar 

We consider a passive scalar T(k, t ) being advected by incompressible turbulence 


^ + Hok 2 ]T(k,t) = -ik a J d 3 ju a (k-j,t)T(j,t) (1) 

with the turbulent velocity field u(k, /) being determined from the Navier-Stokes equation 

[J^ + V 0 k 2 }u a (k,i ) = MaP-y{k) J d 3 juyQ,t)u 1 (k- j,0 + fa(k,t). (2) 

Summation over repeated subscripts is understood, and 

M a p y {k) = k 0 D tn (k) + k r D a p{k) i and D a0 {k) = 8 aP - k a kp/k 2 . (3) 

Here fio is the molecular diffusivity, no the molecular viscosity and f a a random forcing term. 
The forcing correlation is given by 

< >= D 0 k-yD af} (k)8(k + k')S(t-t') (4) 


where D 0 denotes the intensity of the forcing, and y is an appropriately chosen exponent so 
as to recover the Kolmogorov energy scaling in the inertial range (y = 3). 
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2.1 An outline of the recursive RNG procedure 

Since the details of the recursive RNG procedure for Navier-Stokes turbulence have been 
presented before 6,10-11 , we just briefly outline the steps here: 

(i) The subgrid wavenumber region ( k c , kj) is partitioned into fV-shells 

k c = kiv < kpj - 1 < .... < k\ < ko = kd (5) 

where k c is the wavenumber separating the resolvable from the subgrid scales and kd 
is the Kolmogorov dissipation wavenumber. k n — f n k 0 , n — 0,...,N, where / is a 
factor, 0 < / < 1, measuring the coarseness of the subgrid partitioning. The limit 
/ — > 1 corresponding to a differential partitioning of the subgrid region (N — ► oo). 

(ii) The subgrid modes for the first shell, k\ < k < ko , are eliminated from the resolvable 
scale equation by the solution of the subgrid scale equation. 

(iii) A subgrid scale average is performed over the resultant resolvable scale equation. This 
will result not only in the introduction of the subgrid scale energy (or equivalently, 
forcing) spectrum, but it also results in a new triple nonlinearity and nonlocal eddy 
damping function in the resolvable momentum equation ( k < k c ). 

(iv) The above steps are repeated for each successive subgrid shell until all the subgrid 
scales have been removed. 

(v) Since the subgrid scales evolve on a faster time scale than the resolvable scales, a 
multiple time-scale analysis can be performed to simplify the eddy damping function. 
The resultant eddy viscosity is a fixed point of an integro-difference recursion relation. 

(vi) The recursion relation for the eddy viscosity and the renormalized Navier- Stokes equa- 
tion are rescaled. 

It should be emphasized that there are two singular limits 7 : / — > 1 and k —> 0. A careful 
analysis must be done regarding these two limits and the associated averaging operations. 
We will address this issue here in the present paper. 

2.2 An asymmetry in the renormalized passive scalar equation 

The details of the implementation of the recursive RNG procedure to the advection of a 
passive scalar is a straightforward generalization of that for Navier-Stokes turbulence (see 
e.g. Hossain 13 ) and so will not be presented here. Here we comment on a symmetrization 
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procedure that is standard when dealing with Navier-Stokes turbulence, but which can not 
be applied to the passive scalar problem. 

Consider the removal of the first subgrid shell and use the usual notation 


and 



if k\ < k < k { j ; 
if k < k\ . 


( 6 ) 


7Y], _ f t) if < k < k 0 ; 

HK, *)- \T<(k,t) if k < k x . ’ lj 

We find for k < A,q, the resolvable scale passive scalar and Navier-Stokes equation can be 


written as 


[ ~ + 7 ,< (k, i) = -k„J <^j[<(k- j,i) r< (i.0 + “£( k - j.O r< (j,<) 

+^(k-j,()r>(j,i) + ^(k-j,i)r > (j,<)]. (8) 

and 

[^ + I^ok 2 u^{k,t) = f^(k,t) + M ai3l (k) J d 3 i[u|(j,t)u<(k - j,i) 

+2u|(j, t)u<(k - j, t) + uj(j, t)u>( k - j, <)]• ( 9 ) 

The factor 2 in the Navier-Stokes Eq. (9) arises from the symmetry in the j <-» k-j 
interchange. 

We assume isotropy for both the velocity field and passive scalar, so that the subgrid 
velocity - passive scalar correlations are zero (Lesieur 14 ): 

< « > T > >= 0, (10) 

where < . . . > represents averaging over the subgrid scales. 

On applying the recursive RNG procedure to eliminate the subgrid fields in Eqs. (8) 
and (9), the second term on the RHS of Eq. (8) requires special attention. In Ref. 6, the 
frozen flow velocity field is prescribed, so that this term plays no role in the renormalization 
procedure. However, here we are considering a passive scalar field being advected by a 
turbulent velocity field which is itself determined from Navier-Stokes turbulence. As a result, 
this term must be treated on the same footing as the others. For this second term on the 
RIIS of Eq. (8) we have the "wavenumber restrictions, j, k < k i while k\ < |k — j| < 
&o is in the subgrid scale. This can only be achieved when k and j are located near k \ . 
Thus with the wavenumber restriction discussed above is different from 

u aCb t)T < (k - j, t) with k, |k — j| < Aq and k x < j < k 0 . Hence one can not interchange the 
variables j and k-j by standard symmetry arguments as done by Hossain 13 . 
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2.3 The renormalized Navier-Stokes and passive scalar equations 

After removing the n th subgrid shell, the renormalized passive scalar equation takes the form 


j t + fin (k)k 2 ]T(k,t) = -ik a J d 3 ju<(k-j,t)T<(k,t) 
/ d 3 jd 3 j' —^— u <( k - j,o«g(j 


-lk Q t, f d 3 jd 3 j'— 

h=\ J V n- 


M afil (k- j) 


(|k j I) |k — j S 2 e 


7^u<(¥,t)u<(k-j - j , ,f)7 ,< (j,0 


( 11 ) 


where the last term on the RHS differs from that of Hossain 13 . 
The restriction on the wavenumbers are the following: 


kjy < j < k/v-i in the second integral 
k N < |k - j| < &;v-i in the third integral. 

The other wavenumber constraints are as indicated by the superscript on the fields u and T. 
The renormalized Navier-Stokes equation has the following form 11 : 


\d/dt + v n (k)k 2 ]u a (k,t) = /«(k,l) + M a p y (k) J cPju%(j,t)u<(k - j, t) 
+2M„ a ,(k) £ / - j',t)«<(k - j,() (12) 

h-\ J V n -h\J)J 

where j is restricted to the subgrid shell in the second integral. Again, all other wavenumber 
constraints are as indicated by the superscript. 

2.4 Recursion relations for eddy viscosity and diffusivity 

Although the second term on the RHS of Eq. (8) contributes a new triple nonlinear term, 
it does not make any contribution to the renormalized eddy diffusivity in the momentum 
equation in the process of removing the next .subgrid shell. The reason is the following: 

term ~ M a/h { k - j) < u>(j', t)u>(k - j - j',<) > T<{ j,i) ^ 0 (13) 

since the ensemble average will generate a delta function £(k— j) while k— j is in the subgrid 
range. This is impossible, and so this second term can not contribute to the eddy diffusivity. 

After the removal of the (n + l) </l subgrid shell, the spectral eddy viscosity in the renor- 
malized momentum equation is determined by the recursion relation 11 



Vn +l(k) = l' n (k) + 81 / n (k) 


( 14 ) 


where 


and 


8v n (k) = 


Do A f , 3 - I(fc,j, 9 )lk-jr 
k 2 f^ Q J J i'h,(j)j 2 i'(\k - j|)|k- j| 2 


(15) 


L{k,j,q) 


kjj 1 ~ * 2 )N 2 ~ fcj] 

<7 2 


(16) 


with k - j = kjz. This difference equation, after rescaling, has been solved by Zhou et al. 11 
and fixed points were readily determined for finite f < 0.7. However, it was very difficult to 
determine fixed points for finer subgrid partition factor / > 0.7. In Sec. 4 we shall pass to 
the differential subgrid limit / —> 1 and determine an ordinary differential equation (o.d.e) 
for the renormalized eddy viscosity over the entire resolvable scale which can be readily 
integrated. 

In a similar fashion, the spectral eddy diffusivity in the renormalized passive scalar equa- 
tion can be shown to be given after the removal of the (n + l) t/l subgrid shell, by 


^n+l (k) = fin(k) + 8[i n (k) 


(U) 


where 


Sfi n (k) 


kgkp f j3 • 

k 2 hi 


• D a p{k - j)Q(\k — j|) 


Pn-h{j )j 


1 2 


(18) 


The renormalized eddy viscosity and diffusivity are defined as the fixed point of these 
recursion relations. 


2.5 Rescaling of the recursion relation and momentum equations 

From the self-similarity properties of the forcing and energy spectrum in the subgrid range, 
we expect that the viscosity // n + 1 to be simply related to v n for large n, while the diffusivity 
fi n +\ is simply related to fi n . A rescaling can be performed on the recursion relation. In 
particular, consider 


and define 


k — > k n +-\k 


u n {k) = kl y +' )l3 v n (k n+ , k) for k< 1 


(19) 


( 20 ) 
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Unless mentioned otherwise, we drop the tilde notation on the wavenumber, and note 
that 0 < k < 1. The recursion relation for the eddy viscosity becomes 


0 n+1 (k) = f^ /3 [i> n (fk) + 8u n {fk)] (21) 

where 

Do / ,3. I(fc,j,< 7 )[k- j|- y 

k 2 J ^n(/j)j^n(/|k-j|)|k-jP 

" Do /■ 3 h ( y+l y 3 £(fc,j,<7)lk ~ j| V 

tx k 2 J |k-j|)|k-j|U 

and the summation term arises from the triple nonlinearity induced by the recursive RNG 
procedure. The recursion relation for the eddy diffusivity is 


( 22 ) 


Sv„(fk) = 


An +,(*) = /"” + l,/2 [A, I- «An(/*)l 


(23) 


with 


fifinifk) 

n 

+ £ 

/i=i 


_ f- (m+l) [k a k p r D a 0 (k-j)\ k-j|-(-+ 2 ) 

4tt l P J M/»J a 

kgkj f -h(m+\)/2 f ,3 - Dapik - j)|k - j|- (m+2) - 
J J fin-h(f h+l j)j 2 


(24) 


where again the second term on RHS of Eq. (24) arises from the induced triple nonlinearities. 
We have introduced the parameter m = 5/3. These equations are valid for any k in the 
resolvable scales: 0 < k < k c . In the limit A: — > 0, the triple term contribution —> 0, as will 
be shown in Sec. 4.1. 

The final renormalized passive scalar equation is 


+ fi(k)k 2 jr(k,t) = -ih„ J ~ j,i)7' < (k, 4) 

Pur I W^u<(k-i,t)u<Q-j',t)T<(j',t) 

H(k c )k c 2 J 2 

u<(k-j -i',t)T<(j,t) 

v{k c )k c 2 J |k-j| — 

while the final renormalized Navier-Stokes equation is 

[d/d t + v{k)k 2 )u a {k,t) = / a (k,i) + M a0 p k) J d 3 ju^(j,t)u<(k - j ,t) 

+ 2 t±rM a(3l (k) I d 3 jd 3 j' j - j', *)u<(k - j, t). 

v{k c )k c * J J 3 


(25) 


(26) 
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3 Galilean invariance of the renormalized Navier- Stokes 
and passive scalar equations 

Before we proceed further, we turn our attention to the question of the Galilean invariance 
of the renormalized Navier-Stokes and passive scalar equations (25) and (26). The impor- 
tance of Galilean invariance in turbulence modelling has been emphasized by Speziale 12 . 

To be consistent with the basic physics, it is required that the description of the turbu- 
lence be the same in all inertial frames of reference. The appearance of the triple nonlinear 
term, which is a function of the resolvable scales velocity fields, makes the property of the 
Galilean invariance of our recursive RNG procedure not transparent. We now show that 
both the renormalized Navier -Stokes equation and the renormalized passive scalar equation 
are Galilean invariant. 

3.1 Galilean invariance in Navier-Stokes equation: A review 

The Galilean transformation is 


x — > x* — u 0 <* t -> t* 


(27) 


Thus, one has 


u = u* - U 


0 ) 


d_d 
dx dx * 


d_ 

dt 



d 

dx* 


(28) 


While the Galilean transformation for the Navier-Stokes equation in physical space is 


trivial, the Galilean transformation in wavenumber space is less obvious, due to the lack of 
differential operations. For convenience, we first review how Galilean invariance is preserved 
for the Navier-Stokes equation in the wavenumber space. 

Under the Galilean transformation, the LHS of the Navier-Stokes equation [cf. Eq. (2)] 


becomes 


+ U w ikUU k*, i) + j/„fc* 2 [-t/o^(k-) + «k', <)] 
dt* 

= + U ofi ik;u* a (k\t) + v 0 k* 2 u* a {k*,t) (29) 

at* 

where in the last step, we have used the the 6 function property k* 2 S( k*) = 0. 

Under the Galilean transformation, the RHS of the Navier-Stokes equation [c.f. Eq. 2] 
becomes 
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( 30 ) 


AW*') j d 3 j[u' e (r,t) - U w 6 ( j*)][t,J(k' - j',i) - £M(k* - j*)] 

= m„„,(v) J d 3 ru' e (j",t)u;(k" -i‘,t) + iU„k;u’jk’,t) 

where we have used the property of the 8 function, the incompressible condition, and 


Mafo(k*)Uo0U*(k*,t) = U 0 pk*u:(k\t)/(2i) (31) 

Thus, as expected, the Navier-Stokes equation is invariant under a Galilean transforma- 
tion due to the cancellation of the second term on the RHS of Eqs. (29)-(30). 

3.2 The renormalized Navier-Stokes equation under a Galilean 
transformat ion 

To show that the renormalized Navier-Stokes equation is invariant under a Galilean trans- 
formation, we need only consider the recursive RNG induced triple nonlinear term, denoted 
by N St' 


NSt = 2AW*) J -j',tM',t)u y (k~i,t) (32) 


u(k c )j 

It is important to note that j is in the subgrid. 

Under a Galilean transformation, Eq. (32) becomes 


NS* t = 2 M alh (k*) J d 3 j*d 3 i"(^) 4/3 ^^ [vCT - j", t ) - Uop'S(i* - j'*)] 




K(T.*) - ^anK(k*-jv) - - j*)] 


(33) 


Since j* is in the subgrid scale, while j'* and k* are in the supergrid, £(k* — j*) and 
d>(j* — j'*) can never be simultaneously satisfied. As a result, 


NS* t = 2 M at3l {k*) J d 3 j*d 3 j'*( J -) 4/3 - 


j \-1/3 Maf3-y{j ) * /.* •/ 


K,(y*,o-t/o 7 ^(r)K(k*-r,f) 


*( 
7' 


(34) 
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Now only one term in Eq. (34) could violate the Galilean invariance of the renormalized 
Navier-Stokes equation. However, 

/ - r, 0 - t) (35) 

This is not permissible since Upi = u^* and j* is restricted to the subgrid. Thus NSt — 
NSj. Hence the triple term is Galilean invariant. 

3.3 Galilean invariance in the renormalized passive scalar equa- 
tion 

The renormalized passive scalar equation has two triple nonlinear terms. The proof of the 
Galilean invariance of the renormalized passive scalar equation proceeds in a similar manner 
to that for the renormalized Navier-Stokes equation. 

For the first triple nonlinear term, labelled PStu after a Galilean transformation, w T e 
found 


PS? n ~ J d 3 j*d 3 j'*[u<*{ k* 


r,t) + u 0o s(k* 


j*)K* (j' - j“,<) + UoeStf - r))7’ < *(i,() 

/ <e (k- -j'')r < -(j.i)(36) 


since j* is in the subgrid while k*,j r are in the resolvable scale. Thus the £(k* — j*) and 
£(j* — j'*) can never be satisfied simultaneously. 

The second triple nonlinear term has the following structure after the Galilean transfor- 
mation 


ps- T2 ~ J c(>j-d>nu?(r,t) + u w s(r)Hu<‘(k- - j- -r.o + v„y;(k* - y 

- J <fV<f , uJ'Cj'VX(k* -j- -j'-,()r<(j,i)(37) 


where the last two steps follow from the wavenumber constraints. k,j,j' are in the resolv- 
able scales while Jk — j | is in the subgrid scale. Specifically, the first step follows since 
<5(k* — j* — j'*) can never be satisfied. The second step follows since <f>(j ,, ‘) would force 
u^*(k* — j* — j'*) — > u^*(k* — j*). This is not permissible since |k — j | is in the subgrid 
while u*, by definition, is in the resolvable scale. Thus, the renormalized passive scalar 
equation is also Galilean invariant. 
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4 Differential equations for the renormalized eddy 
viscosity and diffusivity 

The differential limit, / — » 1, is singular and has been discussed recently 7 . In particular, it is 
related to the assumption of local versus non-local interactions in k. In this section we will 
calculate the eddy viscosity and diffusivity under the differential equation limit for recursive 
RNG. 

For recursive RNG we will find that the differential equations hold throughout the resolv- 
able wavenumber range 0 < k < k c . This should be contrasted with e - RNG eddy viscosity 
differential equation which is valid only in the k — * 0 limit 8 . 

4.1 The distance interaction approximation, k — » 0 

Consider the resolvable scale Navier-Stokes equation, Eq. (9), 

r\ 

— + v 0 k 2 ]u<(k,t) = M aP ^(k) J d 3 j[u^{j,t)u<(k- j,t) 

+2«|(j,<)«?(k- j,<) + w|(j,0“?(k-j,01 

The first and third terms on the RHS of (9) are symmetric in j and |k — j | in terms of 
their respective wavenumber constraints in wavenumbers. As a result, the distant interaction 
limit k — > 0 has no effect on the existence of these terms which will give rise to the standard 
quadratic nonlinearity and eddy viscosity, respectively. However, the second term on the 
RHS of (9) has the following constraint: j is in the subgrid while |k — j | is in the resolvable 
scales. Specifically, the consistency condition requires that, for small &, j satisfies 

j > k c and j < k c + kz. (38) 

Since \z\ < 1, the range of integration must be O(k). 

Thus, the second term on the RHS of Eq. (9) can not contribute in the limit k — » 0 since 
the integrand is bounded. A similar conclusion can be drawn for the second and third terms 
in the renormalized passive scalar equation. Hence, the triple terms will not contribute to 
the renormalized momentum equations and recursion relation for the transport coefficients 
in the distant interaction limit, k —> 0. However, they will contribute to the renormalized 
Navier-Stokes and passive scalar equations for 0 < k < k c . 
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4.2 Test of the conclusion in §4.1 from numerical simulation 
databases 

The conclusions of the last subsection can be tested directly using numerical simulation 
databases. Indeed, energy transfer and eddy viscosity can be analysed using results from 
numerical simulations by introducing an artificial cut at a wavenumber k c that is smaller 
than the maximum resolved wavenumber k m of the simulation. With this fictitious separation 
between the subgrid and resolvable scales, it is possible to evaluate the effect of the subgrid 
k c < k < k m on the resolvable scales k c . To facilitate comparison with the recursive RNG 
analysis, we consider separately the contribution to the energy transfer and eddy viscosity 
from the second and third term on the RHS of Eq. 9. We form an energy equation from 
the momentum equation and introduce the following notation: is the spectrum of 

energy transfer to mode k resulting from interactions between modes with wavenumbers less 
than k c ; T >< (k ) and T y> (k) represent similar contributions from interactions with one or 
both modes above the cutoff k c , respectively. The equivalent contributions to eddy viscosity 
in the energy equations are u >< (k) = —T >< {k)/2k 2 E(k), and u >> {k ) = —T >> (k)/2k 2 E(k). 

To determine the behavior of the energy transfer and eddy viscosity v >< (k) and u y> (k) 
we measured them in flow fields obtained from numerical simulations on 12b^ meshes of 
forced turbulence. The forced flow dataset was generated by Chasnov lj in a large-eddy 
simulation (LES) of the Kolmogorov inertial range, using a subgrid model derived from the 
stochastic equation that is consistent with Eddy-damped-quasinormal Markovian (EDQNM) 
approximation. 

In Fig.l, we present a numerical measurement of It demonstrates that the second 

term on the RIIS of Eq. (9) does not contribute to the energy transfer process as k — > 0, 
consistent with our analysis. 

4.3 The differential equation limit, / —*■ 1 

We now derive the differential equation for the transport coefficients for finite fc, 0 < k < 1, 
where the wavenumbers are normalized with respect to the cutoff wavenumber k c = 1. The 
o.d.e. in the distant interaction limit ( k — 0) will be derived in the next subsection. After 
the rescaling, we rewrite the recursion relation in the form 

*»+,(*) - / ( ' +,|/ V(/fc) = (39) 

F or f _> i ? the number of iterations n — > oo. Similarity considerations lead to 

u n+ i(k) -> i/(fc), n — > oo. (40) 
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Let A = 1 — /. the LHS of Eq. (39) becomes 


v{k) - [1 - A]< y+l V 3 v[k(l - A)] — ► A 


dv{k) y + 1 


+ 


-v{k) + 0(A)] 


(41) 


dk ' 2 

As noted earlier 7 , the partial average of Rose 6 must be employed in order to insure the 
existence of the differential limit. The partial average is introduced since the distinction 
between the resolvable and subgrid scales become fuzzy in the limit of a differential subgrid 
partitioning, / — * 1. 

Following Rose 6 , we first change the variable from j,z to j,q = |k — j| , with djdz = 
(qfkj)djdq, so that the RHS of Eq. (39) becomes 


+ 


8u(k) = J djdq (~ ) 


L(k,j,q) 


kj t/(j)k 2 u( |k -j|)|k-j| 2 |k-j|v 
L(k,j,q ) j 


l 


dq 


(jLRH-R/s 

kj'v(j)k 2 v (|k - j|)|k - j[ 2 |k — 

L(k ' '•«)_ + A jT 


<?<i+fc k z u 2 {\)q 


<i<i+fc 


Pv 2 ( 1) 


(42) 


where 


L{k , 1 ,q) = fc(l - z 2 )(& - zq 2 )/q 2 (43) 

and 


L(kJ,l) = kj(\ - z 2 )(kj - z). (44) 

Here one has set the coefficient D' Q = 2ttD 0 = 1 (Zhou et al. 11 ). 

As a result, the fixed point renormalized eddy viscosity v(k) is determined from the o.d.e. 
at 0(A) 

= -JfJjlMk) + B„(k)] (45) 

where 


A v {k) 


u 


l+k d kiX 1 ,q) 

nV- 1 


Here, 2 : is evaluated at j = 1 and q = 1, respectively in the L(k,j,q) expression. 
The fixed point o.d.e. for the eddy diffusivity is that given by 


(46) 

(47) 
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(48) 


kdJ ir + 


■K(t) + B,(k) ] 


<g<l+fc 


, ... 1 [ , sin 2 (fc, < 7 , 1 ) 

^(k) “ 2k h< q <\+k q <7 m+1 ’ 

(Jfc) = T^r [ dj sin 2 (fc, l,i)j (m_1)/2 

2k J\<j<i+k 


and sin 2 (fc,j,tf) is the square of the sine of the angle defined by the k and q legs of the 
k,j,q wave- vector triangle. Note that Eqs. (48)-(50) are identical with the o.d.e. of eddy 
diffusivity which Rose 6 derived for a prescribed frozen velocity field. 

4 4 Differential equations in the k — * 0 limit 

In the k -> 0, we have seen that the triple term does not contribute to the eddy viscosity. As 
a result, the recursion relation will now contain only the usual quadratic contribution.. We 
further simply the analysis by taking the standard subgrid linear propagator 11 G h ' (|k - j|) = 

[£ + ^|k-j|)]~G / 7 1 (|j|) asfc-^0. 

The limits of the integration are given by 


1 < j - kz < 1//, 1 + kz < j 

< 1 /f + kz 

(51) 

Thus, the RHS of Eq. (39) becomes 



8u n {k) =M-A-B 


(52) 

where the integral limits for these terms are 



f \/J r\ 

J dj J 1 dz f or 

M 

(53) 

fl r\+kz 

j dz J , d J f or 

A 

(54) 

f° ,. r 

dz dj for 

J - 1 J\/f+kz 

B. 

(55) 


Terms A and B are the corrections to the symmetric term M. They are important for 
a finite bandwidth /. However, it is easy to show that A + B = 0 for / - 1 m the k -» 0 
limit. Hence 


/> -'«•* y - "i - « 


8 gj, 
‘I5i/ 2 (1) 
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while the LHS of Eq. (39) yields 


dk 

since is bounded as k 


u dv{k) y+1 , TN y + 1 /fX 
k - h —z—v(k) -» — v(k), 


0. 


as 


k 


0 


v{k — > 0) = 


3 8 D' 0 

y+ 1 15 // 2 (1) 


Again, one sets the coefficient D' 0 = 2wD 0 = 1. 


( 57 ) 


A similar anlysis can be performed on the fixed point o.d.e for the eddy diffusivity. This 
was not performed by Rose 6 who did not consider the k -> 0 limit carefully. Again, as k -> 0, 
the triple term will not contribute and we find that the corresponding Sft D (k) term has the 
limiting form 


6y D (k) 


1 

1 i 

CM 

[jU/zl 


J dz [ l + (j)y](i -- 2 )p=, 


as k — ► 0. Hence 


A 2^i) />(*-**) 


(58) 


/r(^’ — * 0) 


2 1 
y + i^i)' 


(59) 


4.5 The momentum equation eddy viscosity and diffusivity 

The o.d.e’s, Eqs. (45) and (48), for the momentum equation eddy viscosity and diffusivity are 
readily solved and shown in Fig. 2. We observe that both the eddy viscosity and diffusivity 
have a similar plateau structure as k — > 0. Notice that the eddy diffusivity plateau is not 
obtained in the original numerical calculation of Rose 6 . As k — > k c , eddy viscosity displays a 
weak cusp like behavior while the eddy diffusivity decreases monotonically as k — * k c . These 
curves, as A: > k c , are similar to that of Zhou et aid 1 and Rose 6 . 


5 The turbulent transport coefficient in the second 
moments 

The concept of the spectral eddy viscosity and diffusivity are introduced in the second 
moments 3,4 . Thus, the momentum equation spectral eddy viscosity and diffusivity are only 
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a partial contribution to the total transport coefficients. Indeed from our numerical mea- 
surement in Fig. 1, we expect that the triple nonlinear terms will contribute to the energy 

transfer when k is near k c . . 

We consider the contribution of the triple nonlinear term in the renormalized eddy vis- 
cosity to the eddy viscosity first (Zhou and Vahala' 6 ). The second moment for the velocity 

field is defined as 


U a p( k, f) =< n a (k, t)u«(-k, I) > • 


(60) 


The time evolution of t/ a /?(k,£) 


= -2 v{k)k 2 U a p{k,t) + 2 < f a { k,t)up(-k,t) > + r 5?( k >0 + ( 61 ) 

dt 

In this equation. T^k, t) is the standard energy transfer from the quadratic nonlinearity. 
In contrast, Tj„(k, i) = -2 vrW’m is the energy transfer arising from the RNG induced 
triple nonlinearity. It is readily shown that 


v T (k) = 


L 


*+*■ ... |k-j|-»- ! j» +,/3 

d]dz Kk-1) 


(62) 


2i s{k c )ki m+l/2) k2 

In Fig. 3, we see that v T (k) is the major contributor to the cusp-like behavior of the 
spectral eddy viscosity as k -> 0. For k < 1, it has a backscatter of energy from the subgrid 


and resolvable scales. 

We now define Q(k) as the scalar variance 0(0 =< T(-k,t)T(k,t) >. The dynamic 
equation for Q{k) can be constructed by multiplying Eq. (25) with T{-k,t), followed by an 
averaging operation. Again, a quasi-normal approximation is applied to reduce the fourth 
moment to the product of the second moments. Notice that the last term on the RHS of Eq. 
(25) will not contribute to the spectral equation since < up( j', - j - j , t ) >~ ^( k _ J)’ 

a condition that can not be satisfied since k - j is in the snbgrid scale. 

The dynamical equation for the scalar variance is 


^ + p(fc)fc 2 ]0(M) = £ D + E T (63) 

where £ D is the usual transfer function for the passive scalar. E r (fc) is the additional 
contribution from the triple nonlinear term induced by the recursive RNG procedure 


2 T (k) = - 


fi{k c )k 


2k a J d3 - j)up(j - j') >< T(j')T(-k) > 

= - 2 p T {k)k 2 e{k,t) (64) 


m + l/2 
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where 


H T {k) 


| f ,3 . Dgpjk j)Q(|k j|) 

p(/yfc” +1/! * 2 ■/ > 3 - m/! 


(65) 


and the incompressible condition has been used. It is seen in Fig. 4 that Hr(k) is small when 
k is small. However, as k — ► k c , ^r(k) increases rapidly. 

The solution of uj ’(&) is very similar to that of as k — * k c . They are the major 

contribution to the strong cusp in the eddy viscosity found from the Test Field model 3 and 
EDQNM 4 . Furthermore, Vj(k) also contains the backsc.atter of the energy from the subgrid 
to the resolvable scale. This is a major difference between Hr(k) and v T (k). 

Rose 6 discussed the role of the triple nonlinear terms in physical space. He pointed out 
that it represents the possibility of an exchange of scalar eddies between the resolvable and 
subgrid scales. This effect is an inherent property of measurements made on the passive 
scalar system with instruments which have a spatial resolution limited to an eddy width size 
greater than I (k c . 


6 Spectral eddy viscosity, diffusivity and Prandtl num- 
ber 

The spectral eddy viscosity is simply the sum of the contributions from the momentum 
equation and the triple nonlinear term. The result is presented in Fig. 5. It appears that 
our calulation is in qualitative agreement with the closure theory 3,4 and direct numerical 
measurements 14 18 . In particular, it predicts the correct asymptotic behaviors of the eddy 
viscosity as k — > 0 and k —> k c (Kraichnan 3 ). 

Our spectral eddy diffusivity shows a plateau at k — > 0, in good agreement with the 
EDQNM calculation of Chollet 4 . However, the EDQNM calculation is not unique and de- 
pends on the choice of the parameters. Our diffusivity is in good agreement with EDQNM 
when parameters are chosen according to the direct interaction approximation (DIA) 4 . 

The spectral Prandtl number can be easily determined from our calculated eddy viscosity 
and diffusivity (Fig. 5). It is a function of k and has values ranging from 0.72 ~ 0.92. Note 
that our turbulent Prandtl number at k — » 0 limit is very close to that reported by Yakhot 
and Orszag 8 (0.7179). The values of turbulent Prandtl numbers found experimentally 19 
in the boundary layer are in the range of 0.6 ~ 0.8. Hinze 20 and Tennekes and Lu m ley 21 
pointed out that the transfer of passive scalar and velocity fields may be equally effective. 
Thus the turbulent Prandtl number is about 1. However, Lesieur and Rogallo 18 found that 
their spectral Prandtl number only rose from 0.2 at small k to 0.8 near the cutoff. Lesieur 14 
recently found that the turbulent Prandtl number may be much closer to 1 than that of 
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Lesieur and Rogallo 18 . The EDQNM spectral Prandtl number depends on the choice of two 
EDQNM closure scalar adjustable constants 4 - 14 . In the first case, it remains approximately 
equal to 0.6, even in the vicinity of k c . In the second case, it has a plateau value of 0.33, 
and a cusp close to k c , where it rose to 0.6. 

7 Discussion and conclusions 

In this paper we have applied recursive RNG to the problem of the advection of a passive 
scalar by incompressible turbulence. We have clarified the role of the higher-order RNG 
induced nonlinearities and shown that: (a) The renormalized evolution equations are still 
Galilean invariant (i.e., these higher-order nonlinearities do not destroy the Galilean invari- 
ance of the original equations). This is an important property that needs to be preserved in 
subgrid modeling, especially as one proceeds to more complicated flows and boundaries, (b) 
These higher-order nonlinearities do not contribute to the transport coefficients as k -> 0. 

Now the typical byproduct of the recursive RNG methods is a complicated mtegro- 
difference recursion relation to be solved for the eddy transport coefficients 6 - 10 ". This 
recursion relation is a function of the subgrid shell thickness parameter /. Here, we have 
shown how to pass to the differential subgrid shell thickness limit f 1. In this limit, we 
recover an ordinary differential equation for the eddy coefficients - an equation that is very 
easily solved. 

The o.d.e that is derived in recursive-RNG is fundamentally different from that derived 
by e-RNG techniques. For convenience, we have summarized the Yahkot-Orszag derivation 
of the e-RNG o.d.e in the Appendix. In e-RNG, one is forced into taking the k -> 0 limit 8 9 , 
and the independent variable of the resulting o.d.e is actually the cut-off wavenumber k c . In 
recursive RNG, the independent variable is the resolvable scale wavenumber k, 0 < k < k c , 
with a renormalization transformation that permits k c to be fixed. There is no renormal- 
ization transformation made in the Yakhot-Orszag e-RNG formulation. In the limit k -» 0, 
the eddy transport coefficients from both theories are in very close agreement. This is to 
be expected since the higher-order recursive RNG-induced nonlineanties -> 0 as k — ► 0. 
The slight difference in the eddy coefficients (in the k — > 0) between the two theories can be 
attributed to the treatment of k c : i.e., whether one performs RNG rescaling transformations 
(recursive RNG) or not. The important effect of the triple nonlinearities introduced by the 
recursive RNG procedure are clearly seen in the second moment equations - especially for a 
resolvable wavenumber k — *• k c . 

The spectral eddy viscosity, diffusivity and Prandtl number are determined and we find 
good agreement with both closure theory 3 " 4 and direct numerical simulations 14 " 16 . 
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Appendix: Brief summary of the Yakhot-Orszag’s derivation 
of o.d.e. for the eddy viscosity 

In the e- RNG, a small parameter e is introduced into the exponent of the forcing corre- 
lation function, Eq. (4), 

e = 4 + y- d= 1 + y, (-41) 

where d = 3 is the dimensionality 8 . 

The distance interaction 9 is introduced at the outset 8 , k — > 0. Thus, the triple nonlinear- 
ities do not contribute to the recursion relation, Eq. (14). Furthermore, the eddy viscosity 
is a function of the cutoff wavenumber k c and is obviously independent of the resolvable 
scale wavenumber 0 < k <■ k c . To emphasize this, we shall follow Smith and Reynolds 22 and 
denote their eddy viscosity by i/(k — » 0, k c ). 

By variation of the cutoff wavenumber k c , Yakhot and Orszag 8 derived an o.d.e for 

v(k — > 0 , k c ) 

dk ~ = A$u(k -> 0, k c )\ 2 (k -> 0, k c ) (-42) 

where 


4^3 


A* = 


A (t —* 0, k c ) = 


D 0 


(-43) 


(2*) 3 ’ 30 

where S 3 = 47T 2 is the area of a unit sphere in three-dimensions. These equation should be 
compared with that Eqs. (43)-(45) where the resolvable scale wavenumber k is the variable 
in the recursion relation derived from the recursive RNG, and not the cutoff k c . 

For a given boundary condition v(k — > 0, k 0 ), an analytical solution can be obtained from 
Eqs. (A1-A3). 


v(k -» 0, hf - „(k - 0, t-„) 3 = k; ‘ ‘ 

[z 7 T) S e 

where ko is typically in the order of Kolmogorov dissipation wavenumber. Eq. 
to the familiar inertial range form by taking v(k 0 ,&o) — > Vo for k c — > kj. 
Eq. (A4) takes the form 


(-44) 

(A4) reduces 
In that limit, 


where v(k ■ 
turbulence. 


v{k — ► 0, k c ) ~ 
0, A - o) is neglected since v(k 


2 N 3 Do j 1 /3 r 3 -4 3 

(2tt) 3 -I t 2e 

> 0, k c ) » v(k 


1/3 


(-45) 


0, to) in the fully-developed 
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In order to compare with the inertial range theory of Kraichnan 23 , \akhot and Orszag 
set 

k c = k (46) 

for any k € [0, A,' c ]. From the closure theory analysis 3,4 and direct numerical simulation 
measurements 14-18 , eddy viscosity has rather distinct characteristics in the limits k -* 0 and 
k -> k c . As a result, Eq. (A6) is a rather oversimplified approximation 22 . 
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Fig. 1 Forced eddy viscosity profiles determined from LES databases for the fluid velocity 
at one time instant. Vt 0 t[k) — v >> (k) + where v >> (k) arises from measured LES 

nonlocal subgrid energy transfer, and arises from measured LES local subgrid energy 

transfer. It is important to note that v >< {k) -> 0 as k -> 0 and that u >< (k) arises from the 
w<_ ^ > interaction. Note also the cusp bahvior in u y< [k) as k — » k c . 




Fig. 2 A plot of the momentum eddy viscosity, v(k ), and diffusivity, fi{k), as a function 
of the resolvable scales, 0 < k/k c < 1. These profiles are determined from the o.d.e’s for 
recursive RNG in the limit of differential subgrid shell thickness, / -> 1- 




Fig. 3 The drain eddy viscosity uj{k) arising from the triple nonlinearities in the 
differentia] subgrid shell limit in recursive RNG. r = k c /K p is a parameter in the production- 
type energy spectrum, so that E(k) — ► k 4 as k -> 0. I\ p is a parameter that controls 
the location of the peak in E(k). As r increases, this peak in E(k) moves to smaller k. 
Backscatter of energy from the subgrid scales to the large spatial scales is seen for k/k c < 0.4, 
the region in which vrik) < 0. For r > 1, there is a sharp cusp as k — > k c . 




Fig. 4 The drain eddy diffusivity fir{k) arising from the triple nonlinearities in the 
differential subgrid shell limit of the scalar variance RNG evolution equation. The parameter 
r is as in Fig. 3. Notice that there is now no backscatter of scalar variance, since p, T {k) is 
non-negative for all k. There is a strong cusp as k > k c . 
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